Optimal multiqubit operations for Josephson charge qubits 
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We introduce a method for finding the required control parameters for a quantum computer that 
yields the desired quantum algorithm without invoking elementary gates. We concentrate on the 
Josephson charge-qubit model, but the scenario is readily extended to other physical realizations. 
Our strategy is to numerically find any desired double- or triple-qubit gate. The motivation is the 
need to significantly accelerate quantum algorithms in order to fight decoherence. 
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Quantum computing algorithms are realized through 
unitary operators that result from the temporal evolution 
of the quantum system under consideration. Typically, 
these are achieved with a sequence of universal gates [l| 
which act analogously to the elementary gates of digi- 
tal computers. Quantum computers hold the promise 
of exponential speedup with respect to classical comput- 
ers owing to the massive parallelism arising from the su- 
perposition of quantum bits, qubits; for introductions to 
quantum computing and quantum information process- 
ing, see Refs. |2|. Several different physical implementa- 
tions of quantum computing have been suggested; in par- 
ticular quantum computing with Cooper pairs has been 
proposed p(. 

Superconducting circuits Q feature controlled fabrica- 
tion and scalability Q; their drawback is that the leads 
inevitably couple the qubit to the environment, thereby 
introducing decoherence In a superconductor, the 
number of the Cooper pairs and the phase of the wave- 
function constitute conjugate variables. The majority 
of investigations has focused either on the charge regime 
where the number of Cooper pairs is well defined 7] , or on 
the flux regime where the phase is well defined . Qubits 
utilizing current-driven large Josephson junctions have 
been tested experimentally M. The decisive experimen- 
tal progress reported in Ref. has demonstrated that it 
is possible to realize 10 4 elementary quantum gates with 
Josephson-j unction qubits. Here we consider Josephson 
charge qubits. 

In this Letter we propose a method to construct arbi- 
trary two- or three-qubit quantum gates by solving the 
numerical optimization problem of control parameters for 
a Josephson charge qubit register. We show that it is pos- 
sible to numerically find the required control-parameter 
sequences even for nontrivial three-qubit gates without 
employing elementary gates. Recently, it has been sug- 
gested [llj how to solve a similar problem in the context 
of holonomic quantum computation [l2| , where time does 
not appear as an explicit parameter. In the present con- 
text, the time evolution appears through the Schrodinger 
equation. 



FIG. 1: Instead of implementing the three-qubit quantum 
Fourier transform with the help of elementary gates, we de- 
termine a gate that performs the entire three-qubit operation 
with a single control loop. Note that idle time is avoided. 



The motivation underlying the investigation of this ap- 
proach is the need to overcome effects of decoherence. 
The implementation of a quantum algorithm which is 
composed of elementary gates is rarely optimal in execu- 
tion time since the majority of qubits is most of the time 
inactive, see Fig. [I] The decomposition into elementary 
gates works extremely well with classical digital comput- 
ers. However, in the context of quantum computing the 
number of consecutive operations is strictly limited by 
the short time window set by interactions with the en- 
vironment. It is therefore of prime importance to con- 
centrate on the implementations of quantum algorithms 
[13 0, E3 • We consider the construction of quantum al- 
gorithms out of larger building blocks. Whereas careful 
design and manufacturing can significantly increase the 
decoherence time, our scenario can serve to reduce the 
number of the operations needed. 

The Josephson charge qubit utilizes the number de- 
gree of freedom of a nanoscale Josephson-j unction cir- 
cuit. The states of the qubit correspond to either zero or 
one extra Cooper pair residing on the superconducting 
island, usually denoted by |0) and |1), respectively. The 
Cooper pairs can tunnel coherently to a superconduct- 
ing electrode. The charging energy of the qubit can be 
tuned with the help of an external gate voltage, whereas 
tunneling between the states is controlled with the help 
of an external magnetic flux. 

The explicit single-qubit Hamiltonian for the qubit i is 
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FIG. 2: Schematic illustration of three Josephson charge 
qubits with inductive coupling. The adjustable parameters 
include the gate voltages V* and the enclosed fluxes 



To tackle the challenge posed above, we concentrate on 
finding quantum gates numerically. The structure of the 
Josephson-qubit Hamiltonian is such that it is not im- 
mediately transparent how one would actually construct 
even the basic controlled-NOT gate. We accomplish this 
by considering loops j(t) in the control-parameter space 
spanned by {B x (t)} and {B J z (t)}. Therefore, the function 
j(t) is of the vector form 



7 (*)= [Blif] ... B»{t) Bl(t) ... B»{t)\ 



(3) 



where we have assumed a register of N qubits. The tem- 
poral evolution induces the unitary operator 



U = T exp 
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where the standard notation for Pauli matrices has been 
utilized. Here B l z is a tunable parameter which depends 
on the gate voltage, while B l x can be controlled with the 
help of a flux through the SQUID. Note that setting B\ = 
B x = results in degeneracy. At the degeneracy point, 
there will be no change in the physical state of the system. 
In the case of single-qubit gates, it is easy to see from 
this model that any rotation in SU (2) can be performed 
on the qubits. Note that U(2) is not available since the 
Hamiltonian is traceless. In general, we cannot achieve 
U(2 N ) for N qubits since the Hamiltonian of the entire 
quantum register turns out to be traceless. However, the 
global phase factor is not physical since it corresponds to 
a redefinition of the zero level of energy. 

Qubits can be coupled by connecting them in parallel 
to an inductor, see Fig.0 This scenario has the benefit of 
allowing for a longer decoherence time and that of being 
tunable. The resulting coupling term in the Hamiltonian 
between the qubits i and j is then of the form 



coupling 



-CB x B x cF y 



(2) 



where C is a positive parameter depending on the ca- 
pacitances of the qubits and also on the inductance. It 
follows from Eqs. that one can apply nontrivial 

two-qubit operations by simultaneously turning on the 
SQUIDs of the two qubits, although the a^-term will be 
turned on as well. All the other qubits must have their 
SQUIDs turned off. On the other hand, one-qubit (re- 
operations require that all but one SQUID is turned on. 
Note that in the present context it is actually impossible 
to perform independent operations on any two subsets 
of the quantum register due to the inductive coupling. 
Since one must also take into account the decoherence 
mechanism, it is not practical to let most qubits reside 
at their degeneracy point. The question arises whether 
it would rather prove more efficient to try and find some 
scheme of finding larger quantum operations, instead of 
using elementary gates. 
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where T stands for the time-ordering operator and we 
choose h = 1. The integration is performed along the 
path formed by j(t) where the loop starts at the origin, 
i.e., at the degeneracy point. We will restrict the path 
to a special class of loops, which form polygons in the 
parameter space. Thus the parameters vary in time at 
a piecewise constant speed, and none of the parameters 
is turned on or off instantaneously. We further set the 
time spent in traversing each edge of the polygon equal 
to unity. This limitation could be relaxed, in which case 
the length of each edge in time would be an additional 
free parameter. We also set C = 1 in Eq. (J2J). This can 
be achieved by properly fabricating the inductor, but we 
have every reason to believe that the algorithm will work 
for other choices of C as well. Hence, in order to evaluate 
Eq. (O one only needs to specify the coordinates of the 
vertices of the polygon, which we denote collectively as 
Xy. Numerically, it is easy to evaluate the unitary oper- 
ator in a stable manner by further dividing the loop j(t) 
into tiny intervals that take the time At to traverse. If ji 
denotes all the values of the parameters in the midpoint 
of the i th interval, and m is the number of such intervals, 
then we find to a good approximation 



U Xl * exp(-i#( 7m )At) . . . ex P H#( 7 i) At). (5) 

We now proceed to transform the problem of finding 
the desired unitary operator into an optimization task. 
Namely, any U can be found as the solution of the prob- 
lem of minimizing the error functional 



f(x 7 ) = 



\U-U X ^ 
Here II 



(6) 



over all possible values of X 7 . Here || • \\p is the Frobenius 
trace norm defined as \\A\\p = \/Tt[A^A). The number 
of adjustable vertices of the polygon v is kept fixed from 
the beginning. One needs to have enough vertices to pa- 
rameterize the unitary group SU(2 N ). The dimension of 
this group is 2 2N — 1 and there are 27V parameters for 
each vertex. Thus, we must have 2Nv > 2 2N — 1. We 
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use v = 12 for the three-qubit gates and v = 4 for the 
two-qubit gates. Within this formulation the method of 
finding the desired gates is similar to the recently intro- 
duced method of finding holonomic quantum gates [TH . 
Thus we again expect the minimization landscape to be 
rough and we apply the robust polytope algorithm ^| 
for the minimization. 

We concentrate on finding two- and three-qubit gates, 
since one-qubit gates can be trivially constructed with 
the help of Euler angles. A larger quantum gate could 
be performed by factoring it into two- and three-qubit 
operations, and the implementation for these could be 
found numerically. It seems that quantum operations for 
four, five or more qubits could be found with the same 
method, assuming that sufficient computing resources are 
available. However, even in the case of three-qubit gates 
the optimization task becomes challenging and we need 
to use parallel programming. In the parallel three-qubit 
program, since the function evaluations of f{X 1 ) require 
a major part of the computation, we distribute the work- 
load such that each processor calculates the contribution 
of a single edge of the polygon. In addition, one processor 
handles the minimization routine. 

Let us turn to the results. First, we attempt to con- 
struct a gate equivalent to the controlled-NOT, namely 



U = exp 
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(7) 



The phase factor is needed in order for the gate to belong 
to SU(4). It is already hard to see from the form of the 
Hamiltonian how this gate would be carried out in the 
present setting. Figure illustrates an implementation of 
this gate that has been obtained by minimizing the error 
function in Eq. ©; the error is negligible. This example 
clearly illustrates the potential of our method. 

As a second example, we construct the two-qubit 
Fourier transform. The quantum Fourier transform (see 
e.g., Ref. |2|) is given in the case of two qubits by 




FIG. 3: Control-parameter sequences as functions of time 
that yield the gate in Eq. Q which is equivalent to the 
controlled-NOT. The relative error is on the order of 10 -11 
and 100 discretization points per edge were used. 
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FIG. 4: Control-parameter sequences as functions of time 
that yield the two-qubit Fourier transform in Eq. 0. The 
relative error is on the order of 10 -11 and 100 discretization 
points per edge were used. 



Fo = 



1111 

1 i -1 -i 
1-11-1 
1 -i -1 i 



(8) 



Furthermore, we need to multiply this by exp (i^) in 
order to find a gate that belongs to S77(4). Figure 0] 
shows the resulting loop that has been found with the 
help of the algorithm. In general, the optimization task 
for two-qubit gates can be performed quite easily with 
the help of personal computers. However, finding three- 
qubit gates is already quite time-consuming. It proves 
worth the extra effort to do this, though. 



The three-qubit quantum Fourier transform is 



F, = — 
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where uj = exp(i^). Since detfJ^) = i we must set 
U = exp(-^)F 3 such that U e SU(8). As an ev- 
idence of the success of the three-qubit algorithm, we 
have in Fig. [HI plotted the implementation of the three- 
qubit Fourier transform. We conclude from these three 
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FIG. 5: Control-parameter sequences as functions of time 
that yield the three-qubit quantum Fourier transform (mod- 
ulo a global phase). The relative error is on the order of 10 -5 
and 100 discretization points were used. 

examples that it is possible to find far more powerful op- 
timal implementations of multiqubit quantum gates with 
the help of the minimization scheme. 

To further assess the strength of the technique, we 
compare the number of steps that are required to carry 
out the three-qubit Fourier transform using only two- 
qubit gates with the number of steps required when us- 
ing the full three-qubit implementation of Fig. The 
two-qubit implementation requires effectively four 
gates, see Fig. [I] Since these gates would have to be per- 
formed sequentially, we would need five polygon edges 
per two-qubit operation. This results in 20 edges for 
the whole operation. Using elementary gates would re- 
quire far more edges. Our optimized three-qubit Fourier 
transform, though, only requires 13 edges. Since each 
edge contributes the same amount to the operation time, 
we conclude that our implementation is improved. What 
is more, not all multiqubit gates can be decomposed as 
conveniently as the Fourier transform. For them the gain 
is higher. Thus, increasing the amount of classical com- 
puting resources should yield even better results. 

In conclusion, we have described how to efficiently 
construct two- and three-qubit quantum gates for the 
Josephson charge qubit using numerical optimization. 
An immediate strength of the present scenario is that 
one avoids unnecessary idle time during the logical quan- 
tum operations. Since the loops are traversed at a piece- 
wise constant speed, and no fields are instantaneously 
switched, this method of constructing quantum gates 
should be viable from the experimental point of view as 
well. The effect of finite fall and rise times of pulses on 
the quality of quantum gates has been studied recently 
[l^ . Since we do not use pulses but instead interpolate 
along linear paths in the parameter space, such errors can 



be avoided. It seems reasonable to construct large-scale 
quantum algorithms in multiqubit blocks. This can be 
accomplished by optimizing the gate realization with the 
help of classical computers. 

The authors would like to thank M. Nakahara for 
useful discussions. AON and JJV thank the Research 
Foundation of Helsinki University of Technology and the 
Graduate School in Technical Physics for financial sup- 
port; this work is supported by the Academy of Fin- 
land through a Research Grant in Theoretical Materials 
Physics. We also thank the Center for Scientific Com- 
puting (CSC, Finland) for computing resources. 



* Currently at VTT Information Technology, Microsens- 
ing, POB 1207 02044 VTT, Finland; Electronic address: 
antti.niskanen@vtt.fi 

* Electronic address: juhav@focus.hut.fi 

[1] A. Barenco, C. H. Bennett, R. Cleve, D. P. DiVincenzo, 
N. Margolus, P. Shor, T. Sleator, J. Smolin, and H. We- 
infurter, Phys. Rev. A 52, 3457 (1995). 

[2] J. Gruska, Quantum Computing, McGraw-Hill, New 
York (1999); M. A. Nielsen and I.L. Chuang, Quan- 
tum Computation and Quantum Information, Cambridge 
University Press, Cambridge (2000); A. Galindo and M. 
A. Martin-Delgado, Rev. Mod. Phys. 74, 347 (2002). 

[3] A. Shnirman, G. Schon, and Z. Hermon, Phys. Rev. Lett. 
79, 2371 (1997); D. V. Averin, Solid State Commun. 105, 
659 (1998). 

[4] Yu. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 

Phys. 73, 357 (2001). 
[5] J. Q. You, J. S. Tsai, and F. Nori, Phys. Rev. Lett. 89, 

197902 (2002). 
[6] W. H. Zurek, Rev. Mod. Phys. 75 (in print, 2003). 
[7] Y. Nakamura, Yu. A. Pashkin, and J. S. Tsai, Nature 

398, 786 (1999); Phys. Rev. Lett. 87, 246601 (2001). 
[8] T. P. Orlando, J. E. Mooij, L. Tian, C. H. van der Wal, 

L. Levitov, S. Lloyd, and J. J. Mazo, Phys. Rev. B 60, 

15398 (1999). 

[9] J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 
Phys. Rev. Lett. 89, 117901 (2002); Y. Yu, S. Han, X. 
Chu, S. Chu, and Z. Wang, Science 296, 889 (2002). 
[10] D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, C. 
Urbina, D. Esteve, and M. H. Devoret, Science 296, 886 
(2002). 

[11] A. O. Niskanen, M. Nakahara, and M. M. Salomaa, 
Quantum Information and Computation 2, 560 (2002); 
Phys. Rev. A 67, 012319 (2003). 

[12] P. Zanardi and M. Rasetti, Phys. Lett. A 264, 94 (1999). 

[13] J. Siewert and R. Fazio, Phys. Rev. Lett. 87, 257905 

(2001) . 

[14] J. P. Palao and R. Kosloff, Phys. Rev. Lett. 89, 188301 

(2002) . 

[15] X. Wang, A. S0rensen, and K. M0lmer, Phys. Rev. Lett. 

86, 3907 (2001). 
[16] J. C. Lagarias, J. A. Reeds, M. H. Wright, and P. E. 

Wright, SIAM J. Optim. 9, 112 (1998). 
[17] S. Oh, Phys. Rev. B 65, 144526 (2002). 
[18] We assume that one-qubit operations are embedded into 

two-qubit gates. 



